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We have investigated a closed set of equations for the quark propagator, which has been obtained 
earlier within a new, nonperturbative approach to two-dimensional covariant gauge QCD. It is shown 
that this theory implies quark confinement (the quark propagator has no poles, indeed), as well as 
dynamical breakdown of chiral symmetry (a chiral symmetry preserving solution is forbidden). The 
above-mentioned set of equations can be exactly solved in the chiral limit. We develop an analytical 
formalism, the so-called chiral perturbation theory at the fundamental quark level, which allows 
one to find solution for the quark propagator in powers of the light quark masses. Each correction 
~ ' satisfies the differential equation, which can be formally solved. We develop also an analytical 

f-*^ formalism which allows one to find solution for the quark propagator in the inverse powers of the 

. heavy quark masses. It coincides with the free heavy quark propagator up to terms of order l/mq, 

• where tuq is the heavy quark mass. So this solution automatically possesses the heavy quark flavor 

symmetry up to terms of order I/itiq. At the same time, we have found a general solution for the 
, heavy quark propagator, which by no means can be reduced to the free one. 
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> ' I. INTRODUCTION 

m 

The investigation of two-dimensional (2D) QCD in the context of the Schwinger-Dyson (SD) dynamical equations 
of motion has been initiated by the pioneering paper of 't Hooft |lj. He used the free gluon propagator in the 
light-cone gauge, which is free from ghost complications. He used also the large N c (the number of colors) limit 
C*") ' technique in order to make the perturbation (PT) expansion with respect to 1/N C reasonable. In this case the planar 
diagrams are reduced to quark self-energy and ladder diagrams, which can be summed. The bound-state problem 
r"| ■ within the Bethe-Salpeter (BS) formalism was finally obtained free from the infrared (IR) singularities. The existence 
|' of a discrete spectrum only (no continuum in the spectrum) was demonstrated. Since this pioneering paper 2D QCD 
JL , continues to attract attention (see, for example, review and recent papers [H,Q],0| and references therein). Despite 
its simplistic vacuum structure it remains a rather good laboratory for the modern theory of strong interactions, 
c~| which is four-dimensional (4D) QCD [(|. 

In our previous publications 0, @ we have investigated 2D QCD in the arbitrary covariant gauge for the first 
time. In these works a new, nonperturbative (NP) solution (using neither large N c limit technique explicitly nor a 
weak coupling regime, i.e., ladder approximation) to 2D QCD in the covariant gauge is obtained in the context of 
the above-mentioned SD equations, complemented by the corresponding Slavnov- Taylor (ST) identities. It is well 
known, however, that covariant gauges, in general, are complicated by the ghost contributions. Nevertheless, we have 
shown that ghost degrees of freedom can be considerable within our approach 0. The ghost-quark sector contains 
a very important piece of information on quark degrees of freedom themselves through the corresponding quark ST 
identity. This is just the information which should be self-consistently taken into account. In this way a close set 
of equations has been derived for the quark propagator ■ The main purpose of this Letter is to exactly solve the 
obtained system of equations in the chiral limit and to develop analytical methods of its solution in the general case, 
i.e., for the nonzero current quark masses. Let us emphasize in advance that we have found a general solution for the 
heavy quark propagator, which by no means could be reduced to the free one. All of this will provide the necessary 
basis for future numerical calculations as well. 
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II. QUARK SD EQUATION 

The final system of equations, obtained in Ref. |7J for the quantities in the quark sector, are presented by the quark 
SD equation and the quark ST identity as follows (Euclidean signature): 



iyp,0) = zd f ,S- 1 (p)-S(p)T^p,0)S- 1 (p). (2.1) 

For simplicity, here we remove an overbar from the definitions of the renormalized Green's functions, retaining it 
only for the coupling constant g (which has the dimensions of mass) in order to distinguish it from initial ("bare") 
coupling constant. It contains all known finite numerical factors. T M (p, 0) is obviously the proper quark-gluon 
vertex at zero momentum transfer. The Euclidean version of our parametrization of the quark propagator is as 
follows: iS(p) — pA(p 2 ) — B(p 2 ). It is convenient to introduce the dimensionless variables and functions as A(p 2 ) = 
g~ 2 A(x), B{p 2 ) = g~ 1 B(x), x — p 2 jg 2 ■ Performing further some tedious algebra of the 7 matrices in 2D Euclidean 
space, the system (2.1) can be explicitly reduced to a system of a coupled, nonlinear ordinary differential equations 
of the first order for the A(x) and B{x) quark propagator form factors, namely 



xA 1 = -(1 + x)A - 1 - fh B, 
2BB' = -A 2 +2{fhaA- B)B, (2.2) 

where A = A{x), B = B{x), and the prime denotes the derivative with respect to the Euclidean dimensionless 
momentum variable x. For the dimensionless current quark mass, we introduce the notation fho = mo/g. 
The formal exact solution of the system (2.2) for the dynamically generated quark mass function is 

B 2 (c, fho; x) = exp(— 2x) / exp(2x')P{x') dx', (2-3) 



and c is the constant of integration. Not losing generality it can be fixed as c = p 2 /g 2 , where p 2 is some constant 
momentum squared, and 

v(x) = A 2 (x) + 2A{x)v{x) (2.4) 

with 

v(x) = -rh Q B(x) = xA'(x) + (1 + x)A(x) + 1. (2.5) 
Then the equation determining the A(x) function becomes 

du 2 (r) 

¥1 + 2v 2 (x) = -A 2 {x)ml - 2A(x)v{x)fhl. (2.6) 

dx 



A. Quark confinement 

As was emphasized in Refs.0,0, the important observation is that the formal exact solution (2.3) exhibits the 
algebraic branch point at x = c, which completely excludes a pole — type singularity at any finite point on the real 
axis in the x-complex plane whatever the solution for the A(x) function might be. Thus the solution for the quark 
propagator cannot be presented as the expression having finally a pole-type singularity at any finite point p 2 = —m 2 
(Euclidean signature), i.e., 



(2.7) 
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certainly satisfies thereby the first necessary condition of quark confinement, formulated at the fundamental quark 
level as the absence of a pole-type singularity in the quark propagator :9j . It is well known that such kind of unphysical 
singularity (algebraic branch point at x = c) is due to the inevitable ghost contributions in the covariant gauge QCD. 
However, as was explained in Refs. @>@> it will not cause any problems within our approach in order to calculate 
truly NP quantities, such as quark condensate. The absence of a pole-type singularities in the quark propagator as a 
criterion of confinement at the microscopic level is only first necessary condition. The second sufficient condition of 
this criterion, formulated at the macroscopic (hadron) level, is the existence of a discrete spectrum only (no continuum 
in the spectrum) in the bound-state problem within the corresponding BS formalism 1]. Its discussion is obviously 
beyond the scope of the present Letter. 

B. Dynamical breakdown of chiral symmetry (DBCS) 

From a coupled system of the differential equations (2.2) it is easy to see that this system 
allows a chiral symmetry breaking solution only, 

m = 0, A{x) j= 0, B{x) ^ (2.8) 
and forbids a chiral symmetry preserving solution, 

to = B{x) = 0, A(x) + 0. (2.9) 

Thus any nontrivial solution automatically breaks the 75 invariance of the quark propagator, and therefore certainly 
leads to the spontaneous chiral symmetry breakdown at the fundamental quark level (mo = 0, B(x) 7^ 0, dynamical 
quark mass generation). In all previous investigations a chiral symmetry preserving solution always exists. For 
simplicity, we do not distinguish between B(x) and B{x), calling both dynamically generated quark mass functions. 

A few remarks are in order. A nonzero, dynamically generated quark mass function defined by condition (2.8) is 
the order parameter of DBCS at the fundamental quark level. At the phenomenological level the order parameter of 
DBCS is the nonzero chiral quark condensate determined as < qq >o^ —g Jq° dx Bq{cq,x) within our approach (see 
Ref. and Bq(co,x) is explicitly given below in Eq. (2.11)). In general, it can be formally zero, even when the 
mass function is definitely nonzero. Thus the nonzero, dynamically generated quark mass is a much more appropriate 
condition of DBCS than the quark condensate. One may say that this is the first necessary condition of DBCS, while 
the nonzero chiral quark condensate is the second sufficient one. 

C. The chiral limit 

In the chiral limit (m = 0) the system (2.2) can be solved exactly. The solution for the A(x) function is 

A Q (x) = -aT 1 -^ -exp(-x)}. (2.10) 

It has thus the correct asymptotic properties (see Fig. 1). It is regular at small x and asymptotically approaches 
the free propagator at infinity (x — > 00), which can be formally achieved by the two ways: p 2 — > 00 at fixed g 2 and/or 
by g 2 — > as well. Let us note that the last limit is known as the PT one. For the dynamically generated quark mass 
function B{x) the exact solution is 

/•CO 

B$(c ,x) = exp(-2x) exp(2x')Al(x') dx' , (2.11) 

J X 

where cq = p\jg 2 is an arbitrary constant of integration and p\ is some constant momentum squared for the chiral 
limit case. It is regular at zero. In addition, it also has algebraic branch points at x — cq and at infinity (at fixed Co, 
i.e., when g 2 is fixed). As in the general (nonchiral) case, these unphysical singularities are caused by the inevitable 
ghost contributions in the covariant gauges (for general behavior of this solution see Fig. 2). 

As was mentioned above, A (x) automatically has a correct behavior at infinity (it does not depend on the constant 
of integration, since it was specified in order to get regular at zero solution). In the PT limit (g — > 0), the constant 
of the integration cq and the variable x go to infinity uniformly (cq,x — > 00), so the dynamically generated quark 




FIG. 2: The dynamically generated quark mass function as given by Eq. (2.11). 
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mass function (2.11) identically vanishes in this limit, in accordance with the vanishing current light quark mass in 
the chiral limit. Obviously, we have to keep the constant of integration cq in Eq. (2.11) arbitrary but finite in order 
to obtain a regular at zero point solution. The problem is that if cq = oo, then the solution (2.11) does not exist at 
all at any finite x, in particular at x = 0. This is valid, of course, for the general solution (2.3) as well. 



III. NONZERO CURRENT QUARK MASSES 

Let us formulate and develop now the calculation scheme, which gives the solution of the system (2.2) step by step 
in powers of the light nonzero current quark masses, as well as in the inverse powers of the heavy quark masses. For 
this purpose, it is much more convenient to start from the ground system itself, Eqs. (2.2), rather than to investigate 
the general solution (2.3). For this purpose, let us rewrite the ground system (2.2) as follows: 

xA' + (1 + x)A + 1 = -fn B, 
2BB' + A 2 + 2B 2 = 2m AB. (3.1) 

As was mentioned above, we are interested in the solutions which are regular at zero and asymptotically approach 
the free quark case at infinity. Because of our paramctrization of the quark propagator, its asymptotics have to be 
determined as follows (Euclidean signature): A{x) ^z^oo — 1/{ X + ^o)i B{x) ^ x ^oo — fh /(x + m§), and neglecting 
mo in the denominators for light quarks. The ground system (3.1) is very suitable for numerical calculations. 

A. Light quarks 

Let us now develop the above-mentioned analytical formalism, which makes it possible to find solution of the 
ground system (3.1) step by step in powers of the light (u, d, s) nonzero current quark masses, the so-called chiral 
perturbation theory at the fundamental quark level. For this purpose it is convenient to present the quark propagator 
form factors A and B as follows: 

oo oo 

A(x) = B ( x ) = E ™o5„(x), (3.2) 

n=0 n=0 

where it is formally assumed that m^ l ' d ' s ^ <C 1. Substituting these expansions into the ground system (3.1) and 
omitting some tedious algebra, one obtains 

xA' (x) + (1 +x)A (x) + 1 = 0, 
2B (x)B' (x) +A 2 (x) + 2Bl(x) = 0, (3.3) 



and for n = 1, 2, 3, one obtains 



xA' n (x) + (1 + x)A n {x) = -B„_i(x), 
2P n {x)+M n {x) + 2Q n {x) = 2N n _ 1 (x), (3.4) 



where 



p n(x) = B n _ m (x)B' m (x), M n (x) = A n _ m (x)A m (x), 

m—0 m— 

n n 

Qn{x) = ^ B n _ m (x)B m (x), N n (x) = Y A n _ m (x)B m (x). (3.5) 



m— m—0 



Is is obvious that the system (3.3) describes the ground system (3.1) in the chiral limit (m = 0). As we already 
know, it can be solved exactly (see below as well) . The first nontrivial correction in powers of small m is determined 
by the following system, which follows from Eqs. (3.4) and (3.5), and it is 



xA[ + (1 + x)A 1 = -B Q , {B\B' + BqB[) + AqAi + 2B Q B 1 = A B , 



(3.6) 
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where we omit the dependence on the argument x, for simplicity. In the similar way can be found the system of 
equations to determine terms of order in the solution for the quark propagator and so on. 
Let us present a general solution to the first of Eqs. (3.4), which is 



A n (x) = -x~ x e~ x / dx' e x ' B n -i{x'). 
Jo 



(3.7) 



It is always regular at zero, since all B n (x) are regular as well. The advantage of the developed chiral perturbation 
theory at the fundamental quark level is that each correction in the powers of small current quark masses is determined 
by the corresponding system of equations which can be formally solved exactly. 

Let us write down the system of solutions approximating the light quark propagator up to first correction, i.e., 



A(x) = A (x)+m A 1 (x) + ...., 

B(x) = Bo(x)+fh B 1 (x) + .... (3. 



This system is 



and 



A (x) = -x- 1 (l-e- x ), 4,(0) = -1, (3.9) 

Bl{x) = e- 2x / dx' e 2x A 2 {x'), (3.10) 

J X 

= -x- x e- x f dx' e x ' B (x'), (3.11) 
Jo 

B^x) = e- 2x B^(x) [ X dz e 2z A (z)[B (z) - A 1 (z)]. (3.12) 
In physical applications we also need B 2 (x), so we have 

B 2 (x) = B 2 (x) + 2m B (x)B 1 (x) + ... 

= B 2 (x) + 2rh Q e- 2x f dz e 2z A (z)[B (z) - A 1 (z)} + (3.13) 

J ci 

and the relation between constants of integration Co and c\ remains, in general, arbitrary. However, there exists a 
general restriction, namely B 2 {x) > and it should be real, which may lead to some bounds for the constants of 
integration, while x < cq always remains valid. 



B. Heavy quarks 

For heavy quarks (c, b, t) it makes sense to replace m — > ttlq. In this case it is convenient to find solution for 
heavy quark form factors A and B as follows: 

oo oo 

m 2 Q A(x) = fh Q n A n (x), fh Q B(x) = ^ fh Q n B n (x), (3.14) 

and for heavy quark masses it is formally assumed that Tfig' 6 '*-' 3> 1, i.e., the inverse powers are small. Substituting 
these expansions into the first equation of the ground system (3.1) and omitting some tedious algebra, one finally 
obtains 



B (x) = -1, B 1 (x)=0 7 



(3.15) 
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and 

xA' n (x) + (l + x)A n (x) = -B n+2 (x), n= 0,1,2,3,... (3.16) 

In the same way, by equating terms at equal powers in the inverse of heavy quark masses, from second of the equations 
of the ground system (3.1), one finally obtains 

P {x) + Q {x) - N (x) = 0, Pi(a;) + Qi(a;)-JVi(a:) = 0, (3.17) 

and 

P n+2 (x) + Q n+2 (x)-N n+2 (x) = -^M n (x), n = 0,l,2,3,..., (3.18) 

where P n {z) 1 M n (z), Q n {z), N n (z) are again given by Eqs. (3.5). Solving Eqs. (3.17) and taking into account Eq. 
(3.15), one obtains 

A (x) = B (x) = -1, A 1 (x)=B 1 (x) = 0, (3.19) 
so the final system to be solved further becomes 

xA' n (x) + (1 + x)A n (x) = -B n+2 {x), 

P n +2(x) + Qn+2(x) - N n+2 (x) = -^M n {x), n = 0,l,2,3,... (3.20) 

It is possible to show that all odd terms are simply zero, i.e., A 2n+ \(x) = B 2n+ \{x) = 0, n = 0, 1, 2, 3, .... 
The explicit solutions for a few first nonzero terms are 

A (x) = B (x) = -1. (3.21) 



A 2 {x)=x+^, B 2 {x)=x + 1. (3.22) 

A 4 {x) = -x 2 - ^x - y , B i {x) = -x 2 -\x~\- (3-23) 
Thus our solutions for the heavy quark form factors look like 



.•2 



+ ^ T -^ w + ... + D A {x), (3.24) 



m Q m Q m Q 



where 



D ^ = 2^-^- + -« (3 ' 25) 



and 
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1 OO 

B(x) = —y2fhQ n B n (x) 



— ^^ T -^ T + - + D B (x), (3.26) 



where 



m Q m Q m Q 



1 7x + 3 



D B (x) = rr S --—r S -+... (3.27) 



Summing up, one obtains 



A(x) = — — + D A (x), 

X + TOg 

B(x) = r I^— + D B (x). (3.28) 

X + TOg 

In terms of the Euclidean dimensionless variables, the quark propagator is 

iS(x) = xA(x) - B{x). (3.29) 
Using our solutions, obtained above, it can be written down as follows: 

iSh(x) — iSo(x) + xDa(x) — Db(x), (3.30) 
where iSo(x) is nothing else but the free quark propagator with the substitution fho — * ThQ, i.e., 

iS (x) = - X ~ m 2 . (3.31) 

X + TUq 

Since xDa(x) — Db{x) is of order rhg 3 , then Eq. (3.30), becomes 

iS h (x) =iS (x)+O(m Q 3 ), (3.32) 
i.e., it is reduced to the free quark propagator up to terms of order l/ifiq. 

C. Heavy quarks flavor symmetry 

It is instructive to show explicitly that our solution (3.32) possesses the heavy quark flavor symmetry [ToL fTT) . We 
will show that the quark propagator to leading order (up to terms of order 1/uiq) in the inverse powers of the heavy 
quark mass will not depend on it, i.e., it is a manifestly flavor independent to the leading order of this expansion. For 
this purpose, we must take into account that the argument x, which is the dimensionless momentum of the heavy 
quark, contains itself the heavy quark mass mg. In other words, a standard heavy quark momentum decomposition 
should be used, namely 

= mqv^ + fcp,, (3.33) 

as well as 



x = -f^Xf, = 7 (U (mQU Al + y p ), 



(3.34) 
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where v is the four- velocity with v 2 = — 1 (Euclidean signature). It should be identified with the four- velocity of the 
hadron. The "residual" momentum k is of dynamical origin. In these terms the Euclidean dimensionless dynamical 
momentum variable x = p 2 jg 2 then becomes 



x = —nig — 2niQt — z, (3.35) 

where we denote t — (v ■ y) with y^ = k^/g and z = k 2 /g 2 . 

Substituting these expressions into the Eq. (3.31) and taking into account only the leading order term in the inverse 
powers of tuq, one finally obtains 



iS h {v, y) = iS (v, y) + 0( — ), (3.36) 



where 



1 v - 1 

iS (v,y) = — , (3.37) 

vy 2 

which is exactly the heavy quark propagator [llj. Thus our propagator does not depend on ttiq to leading order in 
the heavy quark mass limit, mo — > oo, i.e., in this limit it possesses the heavy quark flavor symmetry, indeed. 



IV. THE GENERAL SOLUTION FOR HEAVY QUARKS 



It is easy to understand that the chiral perturbation theory at the fundamental quark level developed for light 
quarks in subsection A completely coincides with the general solution (2.3), complemented by Eqs. (2.4), (2.5) and 
(2.6). We use these equations in the chiral limit as input information in the expansions (3.2). However, things are 
not so straightforward in the case of heavy quarks. Developing the chiral perturbation theory in the inverse powers of 
the heavy quark masses in subsection B, we do not use the general solution (2.3), only the system (2.2) itself. In this 
subsection we will show that the general solution (2.3), complemented by Eqs. (2.4), (2.5) and (2.6), for the heavy 
quark mass function possesses much more information than the direct solution of the system (2.2) on account of the 
expansions (3.14) provides at all. 

Starting from the expansion (3.14) for the A(x) function, which contributes into the quark wave function renor- 
malization only, and using exact Eqs. (2.5) and (2.6), it is possible to show explicitly that it is determined by the 
solution (3.28), i.e., it is 



A(x) = — T +D A (x), (4.1) 

x + TOq 

where the D A (x) function is given in Eq. (3.25). Thus in the case of the A(x) function the straightforward solution 
of the initial system (2.2) completely coincides with exact solution, indeed. 

Unfortunately, things are not so simple for the heavy quark mass function B(x), which should be determined with 
the help of the exact solution (2.3), on account of the solution (4.1). Substituting it first into the relation (2.5) and 
then using the relation (2.4) and doing some tedious algebra, one finally obtains 

1 - 2fhn 2m 2 n 

where 



2x 2x 
D A (x) = D 2 A (x) + 2D A (x) - D' A (x) + 2 D A (x) 

(x + ttIq) {x + rriQY 



+ 2xD A {x)D' A (x) - 2 .^ + 2 ^\ d a (x) + 2(1 + x)D\{x). (4.3) 
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So the general solution (2.3) becomes 



B 2 (c 1 m Q ;x) = 



where 



(1 - 2ft yexp,- 2 «)/;^E<|L^ 



2m,Q cxp(— 2x) 



(x' + m^) 3 



dx' + D B (c, mg;x), 



(4.4) 



Db(c, tuq; x) = cxp(— 2x) I exp(2x')DA(x , )dx' . 

J X 



(4.5) 



The first two integrals can be explicitly expressed in terms of the corresponding integral exponential function Ei(x), 
so one has 



B 2 (c,m Q ;x) = (1 - 2m 2 Q )e~ 



(x + m 2 Q ) (c + m 2 Q ) 



— - 2e~ im Q [Ei(2(x + m Q )) - Ei(2(c + m Q ))] 



+ 2fh 2 e- 2x 



n 2c 



,2 c 



2{c + m 2 Q Y (c + fh 2 Q ) 



+ 2e- Am QEi(2(c + m z Q )) 



2x 2x 



(4.6) 



It is convenient to separate the dependence on the constant of integration c, so after some algebra one obtains 



B 2 (c,m Q :x) 



■ 



Q 



{x + m 2 Q ) (x + m 2 Q ) 2 
+ D B (c,m Q ;x) + D B (c,fh Q ;x) 



2e-^ x+m Q ) Ei{2{x + m 2 Q )) 



(4.7) 



where 



D B (c,m Q ;x) 



1 - 2 ™Q c 2(c 

(c + m 2 Q ) 



-x) + 2e- 2{x+rn ^Ei{2(c + m 2 Q )) 



2mf 



+ ■ 



2(c + rh 2 Q ) 2 (c + m 2 Q ) 



= 2(c-x) 



(4- 



Evidently, by no means the exact solution (4.7) can be reduced to the free quark propagator solution (3.28) for the 
B(x) function, except, maybe, of the asymptotical regime (se below). It is regular at zero. 

It identically vanishes in the PT limit g 2 — ► 0, as it should be, since in this case c = x — > oo uniformly (see 
general solution (2.3)). In the heavy quark mass infinite limit (tuq — > oo and g 2 fixed), the quark momentum also 
goes to infinity (i.e., x — > oo see, for example Eq. (3.33)). In this case the constant of integration c remains finite, 
and therefore the composition D B (c,ffiQ;x) also vanishes. Using further the asymptotics of the integral exponential 
function Ei(z) as follows: 



Ei{z)) = e z 



1 1 , 1 ■ 



z z^ 



z — > oo, 



(4.9) 



from Eq. (4.7) one obtains 



B 2 (c 1 m Q ;x) 



™Q 



2 \2 



+ D B (c, m Q ;x). 



(4.10) 
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Choosing negative sign in the square root, one finally obtains 



B(m Q ;x) 



(x + m 2 Q ) 



1 + 



(x + m 2 Q f 



D B (m Q ;x) 



(4.11) 



where Db(™q;i) does not depend on c and its explicit expression is not important here. Thus, only in the uniform 
limit x, ffig — ► oo the heavy quark propagator may become the free one up to the composition Db(itiq; x), similar to 
Eqs. (3.28). 

It is instructive to present explicitly a few first terms of the expansion (3.14) for the B(x) function by substituting 
it directly into the general solution (2.3), on account of the known already expansion for the A(x) function. Omitting 
all the tedious algebra one obtains 



Bl{x) 



1 



e 



2(c-x) 



B 2 {x) 



1 + e 



2(c- 



>(c-2) 



B^(x), 



(4.12) 



and all odd terms are zero. In order to reproduce the free quark propagator case, one has to go again to the limit 
x — > oo and fixed c. Neglecting then the exponentially suppressed terms and choosing the negative sign in front of the 
square root, i.e., B (x) = — 1, one obtains B 2 (x) = x + 1 and so on. Summing up, one gets the free quark propagator 
solution (3.28), indeed. However, even in this case there is a solution with opposite sign, corresponding apparently to 
the free heavy anti-quark propagator. Concluding, let us remind that it is a general feature of nonlinear systems, like 
the initial system (2.2), that the number of independent solutions is not fixed a priori. 



V. CONCLUSIONS 



We have shown that the quark propagator in 2D covariant gauge QCD reveals several desirable and promising 
features, so our conclusions are: 

1) . The quark propagator has no poles, indeed (quark confinement). 

2) . It also implies DBCS, i.e., a chiral symmetry is certainly dynamically (spontaneously) broken for light quarks, 
while a chiral symmetry preserving solution is forbidden. 

3) . The chiral limit physics (i.e., the Goldstone sector) can be exactly evaluated, since we have found exact solution 
for the quark propagator in this case. 

4) . We develop an analytical formalism, the so-called chiral perturbation theory at the fundamental quark level, 
which allows one to find solution for the quark propagator in powers of the light quark masses. Each correction 
satisfies the differential equation, which can be formally solved exactly. 

5) . We develop also an analytical formalism, which allows one to find solution for the quark propagator in the 
inverse powers of the heavy quark masses. It coincides with the free quark propagator up to terms of orderl/mg. So 
this solution automatically possesses the heavy quark flavor symmetry up to terms of order 1/rfiQ. 

6) . We have proved that the exact solution for the A(x) function coincides with the direct solution of the initial 
system (2.2), obtained by using the above-mentioned expansions in the inverse powers of the heavy quark masses. 

7) . At the same time, the exact solution (4.7) for the heavy quark mass function B(x) by no means can be reduced 
to the free heavy quark propagator. So it is not coincided with the solution obtained by the expansions in the inverse 
powers of the heavy quark masses. 

8) . There is no doubt left that by using the expansions in the inverse powers of the heavy quark masses from the 
very beginning, we are loosing some piece of the important information on the structure of the nonlinear system (2.2) 
itself. So its direct solution, obtained by using the straightforward expansions in the inverse powers of the heavy quark 
masses, has a particular character. Evidently, such straightforward solution does not take into account the response 
of the vacuum, which determines the modification of the quark propagator, while the exact solution (4.7) does take 
this response into the consideration. 

9) . The general solution (2.3) does not demonstrate the principal difference in the analytical structure for light and 
heavy quarks propagators. Also at the fundamental quark level the heavy quark mass limit is not Lorentz covariant. 
That is why in the case of heavy quarks we will use the general solution (4.7) rather than solutions (3.30) and (3.36). 
To take into account the vacuum's response is important even for heavy quarks. 

10) . Our approach makes it possible to calculate physical observables from first principles. All results will depend 
only on the renormalized coupling constant (which has the dimensions of mass) and the corresponding constant of 
integration. A physically well-motivated scale-setting scheme is only needed to fix them. 

Our general conclusion is that 2D covariant gauge QCD implies quark confinement and dynamical breakdown 
of chiral symmetry without explicitly involving some extra degrees of freedom. The only dynamical mechanism 
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responsible for them, which can be thought of in 2D covariant gauge QCD, is the direct interaction of massless 
gluons 0)0- This interaction is a main dynamical effect not only in 2D QCD but in 4D QCD as well, i.e., in QCD 
itself. However, to directly generalize the quark confinement mechanism of 2D covariant gauge QCD to 4D QCD is 
impossible. The problem is that in former theory the coupling constant, having the dimensions of mass, plays the 
role of a mass gap, which was introduced and discussed by Jaffe and Witten (JW) in Ref. [l^- I* 1 latter theory 
the coupling constant is dimensionless, and there is none of the characteristic scales in QCD Lagrangian. If QCD 
confines then such a characteristic scale is very likely to exist, and it is not Aqcd, of course, which can be considered 
as responsible for the nontrivial PT dynamics in QCD (scale violation, asymptotic freedom). A possible way how to 
introduce a mass gap responsible for the nontrivial NP dynamics in QCD has been described in Ref. . Its possible 
relation to the above-mentiond JW mass gap is also discussed there. 
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